Electric field control of the magnetic chiralities in ferroaxial multiferroic RbFe(Mo04)2 
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The coupling of magnetic chiralities to the ferroelectric polarisation in multiferroic RbFe(Mo04)2 
is investigated by neutron spherical polarimetry. Because of the axiality of the crystal structure 
below Tc = 190 K, helicity and triangular chirality are symmetric-exchange coupled, explaining the 
onset of the ferroelectricity in this proper-screw magnetic structure — a mechanism that can be 
generalised to other systems with "ferroaxial" distortions in the crystal structure. With an applied 
electric field we demonstrate control of the chiralities in both structural domains simultaneously. 



Multiferroic materials, in which ferroelectricity and 
magnetic order coexist, are attracting conspicuous inter- 
est as candidates for novel applications in digital storage 
devices [1 . The requirement of having a strong magneto- 
electric effect has focussed the research on compounds in 
which ferroelectricity appears as a consequence of mag- 
netic ordering (so called 'type-II multiferroics'), lead- 
ing to the realisation that cycloidal multiferroics, such 
as TbMnOs j2j and Ni3V208 |3], exhibit an exception- 
ally strong cross-coupling between the different types of 
order. In these compounds the atomic spins rotate 
within a plane that contains the propagation direction 
of the incommensurate modulation ri,i+i which connects 
neighbouring atoms and S^+i. 

A phenomenological theory based on symmetry anal- 
ysis of the magnetic and ferroelectric order parameters 
established that the magnetoelectric coupling has a tri- 
linear formP-[M(V-M)-(M-V)M] in the free energy [4], 
leading to a polarisation given by P cx ri,i+i x (S^ x S^+i) 
[5] . A number of microscopic theories have been proposed 
that respect these symmetry constraints for the canonical 
multiferroic materials [6^, . 

An intense experimental effort and a refined under- 
standing of the symmetry requirements for spin-driven 
ferroelectricity have expanded the range of candidate 
multiferroic materials. An interesting line of research has 
developed specifically on magnetochiral or proper-screw 
systems, in which the atomic spins rotate perpendicu- 
larly to the propagation direction of the screw [8 ; the 
model in ref. 5 predicts no polarisation in this case. One 
mechanism that can lead to electrical polarisation in a 
proper-screw magnetic structure is the coupling to struc- 
tural axiality (ferroaxial coupling) [9 . A crystal structure 
can be considered 'axial' if there exists a structural dis- 
tortion that is unchanged by inversion and makes the two 
senses of rotation (clockwise or counterclockwise) about 
an axial vector A distinguishable. As pointed out in ref. 
[9l the magnetic helicity a can be coupled phenomenologi- 
cally to the electric polarisation P and structural axiality 
A to create a trilinear form, a A • P, which is invariant 



under spacial inversion and time reversal. Very recently, 
"giant" magnetically-induced ferroelectricity has been re- 
ported in the ferroaxial system CaMn70i2 [10 . 

RbFe(Mo04)2 (RFMO, ferroaxial below 190 K) is an 
extremely interesting system to test the interplay be- 
tween ferroaxiality, magnetism and ferroelectricity. Be- 
low Tn ~ 4 K, it orders magnetically in a complex struc- 
ture that possesses both helicity and triangular chirality 
(see below), becoming ferroelectric at the same tempera- 
ture [11 . Yet, its layered crystal structure and exchange 
pathways are sufficiently simple for it to be considered 
the "hydrogen atom" of ferroaxial multiferroics. In pre- 
vious work, it has been suggested that the 120° magnetic 
structure in each layer is in itself sufficient to break in- 
version [11 . Here, we demonstrate that the axiality and 
magnetic helicity also play a crucial role in the onset of 
ferroelectricity. We show that the helicity and triangular 
chirality are in fact coupled together in the free energy 
by the axial distortion, and we present data that show 
switching of both parameters simultaneously with an ap- 
plied electric field. 

RFMO undergoes a structural transition at Tc = 190 
K in which the M0O4 tetrahedra rotate (Fig. [T]), lowering 
the symmetry from P3ml to P3 (ferroaxial point group 
3) [T?, "T3^. Below Tn the Fe spins (one per unit cell) or- 
der magnetically in the ah plane in a 120° structure, and 
rotate in a helix between one plane and the next with an 
incommensurate propagation along c* {qz ^ 0.44, here 
always chosen to be positive). Fig. [TJa) shows a unit 
cell of RFMO with the location of the M0O4 tetrahe- 
dra indicated (oxygen atoms mediate the relevant ex- 
change paths — see below). One can see that there is 
a M0O4 tetrahedron above (or below) each magnetic tri- 
angle. We define triangles as "positive" ("negative") if 
they are associated with a tetrahedron pointing along 
the positive (negative) c* direction [this is indicated by 
the two unit vectors v+ and v~ in fig. [ija)]. In plan 
view, one can rotate the structure in such a way that all 
the "positive" triangles point "up" the page [this orien- 
tation is depicted in Figs III a)-[llc)]. We define a stag- 
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A = A^^O 



FIG. 1: (Color online) The structure and exchange paths of 
RFMO. (a) 3D view showing in-plane (Ji) and vertical (J2) 
paths, as well as the two diagonal paths (Ja and Jb) which 
depend on the ferroaxial distortion. The directions of the 
M0O4 tetrahedra (v^ or v~) are indicated, (b) Shows the 
effect of applying the 3- fold symmetry on the interactions. 
(c)-(e) Illustrate the difference between the Ja and Jb paths 
arising from the distortion parametrized by the axial vector 
A. (Arrows point from 2; = to 2; = 1.) 



gered triangular chirality — a parity odd quantity — 
as at = (Si X S2 + S2 X S3 + S3 X Si) • v/S'^. In other 
words, (Jt = 1 if the spins rotate counterclockwise as one 
circumscribes counterclockwise a "positive" triangle of Fe 
spins (with a v+ tetrahedron at its centre); at = — 1 
otherwise. The magnetic helicity is defined as usual as 
cfh = (S^^o X S^^i) • tqi/S'^ = ±1 for a right- and left- 
handed magnetic screw, respectively. 

There are two possible magnetic propagation vec- 
tors: qi = (1/3, 1/3, (7;^) for atah = 1 and q2 = 
(—1/3,-1/3,^2) for <Jt(Th = —1, leading to distinct sets 
of peaks in the magnetic scattering for the two different 
combinations of chiralities. 

Fig. [TJa) shows a unit cell of RFMO with the loca- 
tion of four exchange paths (mediated by the M0O4 oxy- 
gens) indicated. The super-superexchange terms across 
the prismatic faces, Ja and Jb, are equal by symmetry in 
P3ml but become distinct in P3, so that the difference 
(Ja — Jb) is proportional to the amplitude and sign of 
the ferroaxial distortion [Figsjljc) and[TJd)]. Figs. [TJb) 
andjlje) show the effect of the ferroaxial rotation on the 
exchange paths. Assuming the 120° structure as given, 
the symmetric-exchange energy per Fe ion is 

E = Eo^ ahat( Ja - Jb)^' sin(2^^,) (1) 

where £^0 = — 3Ji6'^/2 + J2S'^ cos{27rqz) doesn't depend 
on the magnetic chiralities. Minimising with respect to 
Qz results in tdin{27Tqz) oc 5'^(Ja - Jb)/(2J2 - {Ja + Jb))- 
We arrive here at our first important result: in the pres- 



ence of a ferroaxial crystal structure, adjacent 120° tri- 
angular magnetic planes will rotate with respect to each 
other, forming a helix, without the need for antisymmet- 
ric or inter planar next-nearest neighbour interactions. 
In a real crystal, both ferroaxial domains will be present 
with roughly equal populations (one with Ja > Jb, the 
other with Jb > Ja). Therefore, the lowest energy mag- 
netic configurations will depend on the 'axiality' of the 
domain in question. Notating the chiralities of a partic- 
ular magnetic structure as (crt,crh) = (±,±), eqn. ([T]) 
implies that {(+,+),( — ,—)} states will be lower in en- 
ergy than {(+,—),(—,+)} if Ja > Jb, or vice versa for 
Ja < Jb- One important consequence is that each of the 
two structural axial domains will order with a distinct 
magnetic propagation vector, say qi for positive axial- 
ity and q2 for negative axiality (the exact combination 
depends on the sign of the magneto-elastic interaction). 
Therefore, each set of distinct magnetic peaks (qi or q2 
peaks) probes a single axial domain. 

Spherical neutron polarimetry is an ideal technique to 
study these magnetic structures and their relationship 
with the electrical polarisation, since it has the ability 
to distinguish domains with different chiralities. In the 
present case the magnetic peaks do not overlap with the 
nuclear peaks in reciprocal space, so it is possible to 
calculate the polarisation of the scattered beam purely 
from the magnetic structure factor M [14]. For this we 
work in the Blume reference frame, in which the X-axis 
is along the scattering vector Q, the Z-axis is vertical, 
and the F-axis completes the right-handed set (see in- 
set to Fig. |2|. Magnetic neutron diffraction is only 
sensitive to the component of M perpendicular to Q, 
which we write = {0^ M±y^ M±z)' For a structure 
with one magnetic ion per unit cell, the Ith moment (in 
the unit cell of the lattice vector R^) is given by = 
/i^o(u T ^v) exp(iqi • R/) + c.c. for (+, +) or (— , — ) struc- 
tures; and is given by /l^^ = /i^o(u±zv) exp(zq2 • R^) +c.c. 
for (+, — ) or (—,+), where u and v are orthogonal unit 
vectors in the plane of the spins. Thus the magnetic 
structure factor M(Q) = pf^^^ {Q){u ± iv) where p 
is a constant, f^^^ is the magnetic form factor, and 
the ± is determined by the magnetic chiralities (see ta- 
ble |l|. Given a fully polarised incident beam directed 
along i, the polarisation measured along j is written Pij 
{i^j = X^Y^Z). For our geometry, the equations reduce 

to Pxx — I5 Pxy — Pxz — J^zy — -^yz — 

Pyy = MljMl - 1; (2) 
Pzz = MlJMl - 1; (3) 
Pyx = P^x = 2^{M^yMl,)/Ml (4) 

where Ml = M± ■ and . = M^iMl-. Thus, 
the sign of the Pyx and Pzx elements, together with the 
position of the satellites (±qi^2) uniquely determines the 
domain population for each of the (at, ah). 

RFMO single crystals were grown by fiux technique us- 
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FIG. 2: (Color online) c-axis polarisation (obtained by 
integrating the pyroelectic current) and intensity of the 
(1/3, 1/3, Qz) magnetic reflection (measured without polarisa- 
tion analysis) as a function of temperature. Inset: the experi- 
mental geometry showing the crystal orientation, direction of 
applied E-field, and scattering vector Q = — k. 



ing high purity (> 99.9%) RbsCOs, FesOs, and M0O3, 
according to the recipe described in JTT . Single crys- 
tal thin platelets (up to 1 cm in diameter) were sepa- 
rated from the flux, their quality and orientation were 
checked using an Agilent Technologies SuperNova labora- 
tory diffractometer, and pyroelectric currents were mea- 
sured and integrated to give the ferroelectic polarisation 
(along c) as a function of temperature (see Fig. |2|. Gold 
contacts were evaporated onto the (0, 0, ±1) surfaces and 
the OrientExpress neutron back-reflection Laue diffrac- 
tometer [15 at the Institut Laue-Langevin (ILL), Greno- 
ble, France was used to mount the sample with a* + b* 
parallel to the F-axis (see inset to Fig. |2|. This allowed 
us to access peaks of the form {hhl). The contacts were 
connected to gold wires with silver epoxy to allow high 
voltage to be applied. An ILL 'orange' cryostat provided 
cooling during the neutron scattering experiment (using 
a fixed neutron wavelength of 0.825 A), which was car- 
ried out using the CryoPad [16 set-up on beamline D3 
at the ILL. 

Initially the sample was cooled below Tc, fixing the 
population of axial structural domains for the rest of the 
experiment. Fig. [2] shows the temperature dependence of 
the intensity of the (1/3, 1/3, 0.44) magnetic peak (corre- 
sponding to the propagation vector qi), which fits the py- 
roelectric data well and confirms the simultaneous onset 
of magnetic ordering and ferroelectricity at Tn ~ 4K. The 
same sample was cooled in an applied electric field of both 
± 7.5 kV/cm and zero field, and the polarimetry compo- 
nents measured. The results for the (—1/3,-1/3,0.56) 
peak are shown in Fig. [sj^a), together with the calcula- 
tions (which contain no free parameters) from equations 



TABLE I: The eight possible contributions of the magnetic 
structure to the scattered intensity, at : the triangular chiral- 
ity; ah' the helical chirality; Q: the position of the peak 
in reciprocal space; G: a reciprocal lattice vector; qi = 
(1/3,1/3,^2); q2 = ( — 1/3,-1/3,^2). The magnetic struc- 
ture factor M is given in terms of the orthonormal vectors u 
and V. 
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([2| to Q. It is clear that the components Pyx and P^x 
couple to the electric field, since in the case of zero field 
cooling these components are both zero. 

Table [l] lists the magnetic peaks and associated mag- 
netic structure factors originating from different chiral 
structures. Peaks with propagation vectors qi and q2 
are present with similar intensities, yielding four satel- 
lites round each reciprocal lattice node; this shows that 
the crystal contains a roughly equal population of both 
axial domains. We examined the behaviour of several 
magnetic peaks, including at least one from each combi- 
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FIG. 3: (Color online) (a) Polarimetry components following 
negative, positive, and zero field cooling (with field strength 
7.5 kV/cm). Solid bars indicate observed values and rectan- 
gles show the calculations, (b) Hysteresis loops in Pzx as a 
function of applied electric field for two magnetic reflections. 
The chiralities, (at, ah) = (i? i) are indicated. 
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nation of chiralities, in each of the field coolings. In zero 
field cooling, the off-diagonal components Pyx and Pzx 
for all peaks are zero. This is only possible if the frac- 
tions of (+, +) and (— , — ) domains (which have equal 
and opposite ^{M}) are equal, i.e. /+,+ = /-,- = 0.5, 
and likewise /+,- = /-,+ = 0.5. For positive field cool- 
ing, the non-zero off-diagonal terms require that only 
two of the chirality combinations be present, so that 
/+,+ = /+,_ = 1 and /_,+ = /_,_ = [Fig. ga)]. For 
negative field cooling, the other two domains are popu- 
lated, so that = = and = = 1 [Fig. 
[4]^b)]. Fig. [sj^b) shows the hysteresis in P^x as a func- 
tion of applied electric field for two magnetic peaks. The 
(1/3,1/3,0.44) peak (resulting from one axial domain) 
switches between (+, +) and (— , — ) at ±V whereas the 
(—1/3,-1/3,0.44) peak (from the other axial domain) 
switches between (+, — ) and (—,+). This confirms that 
we have direct control over the magnetic structures in 
both axial domains by using an applied E-field. 

Since the direction of the electrical polarisation P = 
(0, 0, P^), as deduced by the direction of E, is coupled to 
the magnetic structure, from symmetry considerations 
the magnetoelectric free energy can be written as 



S'^Pz{ciat^C2Aai,) 



(5) 



where the axial vector A = (0,0, A) and ci^2 are con- 
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stants [17]. The two terms in eqn. ([5| are made pro- 
portional to each other by eqn. ([T]) (since the sign of A 
determines the difference in Ja and Jb) which imposes 
that (Jt = — ^cTh- Our results show that the overall tri- 
angular chirality of the entire field-cooled sample is uni- 
form and switches with the electric field, as conjectured 
by Kenzelmann et al. [TT], in spite of the axial domain 
structure. However, as explained above, reversing the 
electric field actually switches both <Jh and <Jt simultane- 
ously within a single structural domain (Fig. [4|. The two 
chiral coupling terms to at and A(7\^ are both of relativis- 
tic origin, since, as we have already shown, domains with 
opposite polarisation have the same symmetric-exchange 
energy. However, the two terms rely on different micro- 
scopic mechanisms; since the second term is proportional 
to the axial rotation, a systematic study of isostructural 
compounds, through measurements or first-principle cal- 
culations, could reveal which of the two provides the dom- 
inant contribution to the development of the electrical 
polarisation. 

We have demonstrated that the axial distortion of the 
crystal structure plays a crucial role in stabilising the he- 
lical magnetic structure of RbFe(Mo04)2, linking trian- 
gular chirality with helicity. By means of neutron spheri- 
cal polarimetery, we determined uniquely the populations 
of domains with each of the combinations of helical and 
triangular chiralities in the two axial domains. By ap- 
plying an external electric field, the domain population 
switches between the pairs of magnetic structures that 
are energetically preferred in each axial domain. The cou- 
pling between magnetic structure and electrical polarisa- 
tion is of relativistic origin, and involves both triangular 
chirality and helicity, the latter term being proportional 
to the axial distortion (ferroaxial coupling). 
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FIG. 4: (Color online) Magnetic structures present in (a) pos- 
itive, and (b) negative field cooling, each of which has two 
contributions (from each axial domain). The direction of the 
ferroelectric polarisation (P) is shown, and the axial distor- 
tion is indicted by the direction of the circular arrows. The 
direction of the M0O4 tetrahedron associated with each tri- 
angle is shown by the ± signs. 
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